Enhanced diffusion and ordering of self-propelled rods 
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Starting from a minimal physical model of self propelled hard rods on a substrate in two di- 
mensions, we derive a modified Smoluchowski equation for the system. Self -propulsion enhances 
longitudinal diffusion and modifies the mean field excluded volume interaction. From the Smolu- 
chowski equation we obtain hydrodynamic equations for rod concentration, polarization and nematic 
order parameter. New results at large scales are a lowering of the density of the isotropic-nematic 
transition and a strong enhancement of boundary effects in confined self-propelled systems. 
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Self propelled particles consume energy from internal 
or external sources and dissipate it by actively mov- 
ing through the medium that they inhabit. Assemblies 
of interacting self-propelled particles (SPP) exhibit rich 
collective behavior, such as nonequilibrium phase tran- 
sitions between disordered and ordered (possibly mov- 
ing) states and novel long-range correlations. Biologi- 
cally relevant systems that belong to this class include 
fish schools, bird flocks [H, bacterial colonies 0| and cell 
extracts of cytoskeletal filaments and associated motor 
proteins (jj. A non-living realization may be a vibrated 
monolayer of granular rods [3] . Collections of SPP have 
been the focus of extensive experimental 0,0, 0] and the- 
oretical studies in recent years. A number of distinct the- 
oretical approaches have proved fruitful for understand- 
ing the complex dynamics of these nonequilibrium sys- 
tems. These include numerical studies of simple mod- 
els inspired by the seminal work of Vicsek 
and phenomenological continuum theories based on 
general symmetry arguments [ll|]. Recent work on de- 
riving the hydrodynamic equations from specific micro- 
scopic models has led to some insight into the orig in of the 
collective behavior of these systems [H, Q3, EJEl, 
An important open question that we address here is the 
interplay between self propulsion and steric effects aris- 
ing from the shape of the particle in controlling the large 
scale physics. 

In this paper we consider a physical model of self- 
propelled hard rods that interact with each other solely 
through excluded volume. The rods move on a passive 
substrate. Self-propulsion is modeled as a nonequilib- 
rium velocity vq along the direction of the rods' long axes. 
The goal of our work is to understand how self-propulsion 
modifies the diffusion processes and the mean-field On- 
sager excluded volume interaction [jjj. Using the tools of 
nonequilibrium statistical mechanics we derive a modified 
Smoluchowski equation that differs from the familiar ver- 
sion for thermal hard rods [ijj in three respects. The first 
and obvious modification is a convective mass flux at the 
self-propulsion speed vq along the direction of orientation 
of the rod. Secondly, self-propulsion enhances the longi- 



FIG. 1: A cap-to-side collision of two self-propelled hard rods 
(the width of the rod is exaggerated for clarity), k is a unit 
vector from rod 2 to rod 1 normal to the point of contact. 
Points on the side of the rods are identified by vectors £i. 



tudinal diffusion constant Dm of the rods, according to 
£)|| — > D\\(l+VQ/kBT). This enhancement arises because 
self-propelled particles perform a persistent random walk, 
as recently pointed out by other authors [HI • Finally, the 
momentum exchanged by two rods upon collision is ren- 
dered highly anisotropic by self-propulsion thus modify- 
ing the Onsager form of the excluded volume interaction. 
This leads to novel anisotropic forces and torques from 
steric repulsion in the Smoluchowski equation. 

These modifications of the Smoluchowski equation 
have dramatic consequences for the properties of the sys- 
tem on hydrodynamic scales. This is illustrated by two 
examples. First, we show that the additional momentum 
transfer from self-propulsion lowers the density of the 
isotropic-nematic transition, thereby providing a micro- 
scopic identification for the physical mechanism respon- 
sible for the enhancement of orientational order observed 
in numerical simulations of motility assays jl| . Secondly, 
we demonstrate that self-propulsion greatly enhances the 
effect of confinement and the role of boundaries. 
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The microscopic model. We consider quasi two- 
dimensional hard rods of length I and thickness 2R con- 
fined to a plane, as shown in Fig. [TJ The i-th rod is 
characterized by the position Yi of its center of mass and 
a unit vector Uj = (cos 6i , sin 6% ) directed along its long 
axis. Each rod free-streams on the substrate, until it 
collides with another rod. The collision results in instan- 
taneous linear and angular momentum transfer such that 
the total energy, linear and angular momenta of the two 
rods are conserved. The microdynamics of the system is 
governed by coupled Langevin equations, 

<9v 



= - T (i, j) Vi + Fui - Ci ■ v, + m (*) , (1) 
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t =-^T(i,j> i -CV + r?f(t) , (2) 



dt 



where Vj = dt^i and u>i — dtOi are the center of mass 
and angular velocities, is the friction tensor, with 
C/3 = C\\UiaUii3 + (±{S a i3 - Ui a Uiff), (r is the rota- 
tional friction, and the mass of the rods has been set 
to one. The second term on the right hand side of 
Eq. |l| describes self propulsion as a center of mass 
force F acting along the long axis of each rod. This 
force is nonequilibrium in origin and arises from an 
internal or external propulsion mechanism. The ran- 
dom forces r/i and rjf describe Markovian white noise 
with correlations (rji a (i) rjjp (£')) = A l a aSij6 (t — t') and 
(vf (*) V? (*')) = A R S i:j S(t-t'). For simplicity we as- 



sume the equilibrium-like form A^, 



and 



A R = 2k B T a ( R /I, with / = £ 2 /12 the moment of iner- 
tia of the rod and T a an effective temperature defined by 
these relationships. Finally, the collision operator T 
generates the instantaneous momentum transfer between 
rods at contact and is given by 

T(l,2) = J J | Vi 2 • k| 6 (-V 12 ■ k) 

xS(T cont )(b 12 -l) , (3) 

where k is the unit normal at the point of contact of 
the two rods directed from rod 2 to rod 1, as shown 
in Fig. [TJ The function T cont (ri, r 2 , £i, £2) is nonzero 
when two rods are at contact and zero otherwise. Here 
£i is a vector from the center of mass^ of the i-th rod 
to the point of contact, & = Siiii ± i?k, where —1/2 < 
Si < £/2 parametrizes the distance of points along the 
axis of each rod from the center of mass and J s ... = 

vi - v 2 + oji x fa - (jj 2 x £ 2 



/*^„ ...dsi. Also, V12 



is the relative velocity of the two rods at the point of 
contact. Finally, the operator bi 2 replaces precollisional 
velocities with their postcollisional values, as obtained 
by requiring energy and momentum conservation. The 
explicit calculation of the T operator is given in [3| . 



Modified Smoluchowski equation. We are interested 
here in the overdamped limit, when inertial effects are 
negligible and the low density dynamics is described by 
a Smoluchowski equation for the the probability distribu- 
tion c (x, t), with x = (r, 6), of rods at a point r oriented 
in the direction 9. The derivation of the Smoluchowski 
equation for self-propelled hard rods can be carried out 
following closely that of thermal hard rods and is given 
in [jjl]. Here, we outline the key steps involved. 

1. First, the noise averaged statistical mechanics of a 
system described by a set of coupled Langevin equations 
is given in terms of the Liouville-Fokker-Planck equa- 
tion governing the dynamics of an N particle distribution 
function [13]. This can in turn be converted into a hierar- 
chy of equations for reduced distribution functions anal- 
ogous to the BBGKY hierarchy for Hamiltonian systems. 
At low density, neglecting two particle correlations, the 
first equation of the hierarchy gives a closed Boltzamnn- 
Fokker-Planck equation for the one particle distribution 
function f(x,p,t), with p = (v, 

2. The probability distribution is c (x, t) = J f(x,p,t). 
In the regime of large friction, the velocities of the 
rods decay to a stationary value on microscopic time 
scales. We use an approximate solution of the nonin- 
teracting Fokker-Planck equation valid in the large fric- 
tion regime, f (x,p,t) = c{x,t) f M (p\6), with f M ~ 

cxp 2 fcir ( v ~ l 'o2) 2 — 2klr -^- j2 ) a Maxwellian dis- 
tribution centered at the self-propulsion velocity vqU. 
With this ansatz, the Bolztmann-Fokker-Planck equation 
can be transformed to a closed equation for the spatial 
probability distribution, c. 

3. To obtain this closed equation we need to evaluate 
the mean force and torque on a given rod due to all other 
rods in the fluid, namely (T(l, 2)vi)m and (T(l, 2)u>i)m, 
where (...) M = / ■■■fM{pi\0i)f M {p 2 \6 2 ). In the ab- 
sence of self propulsion, this average can be readily car- 
ried out and yields the Onsager excluded volume inter- 
action. For finite self propulsion, Jm depends on the 
angular coordinate and hence averaging over velocities 
induces orientational correlations that cannot be incorpo- 
rated exactly. To make progress, we let = Vj — UiVo in 
the calculation of the velocity averages and then neglect 
the coupling between velocity and angular correlations 
by approximating (T(l,2)vi) M - (T(l, 2)vi) M |„ o=0 + 
(T(l, 2)vi)„ , where the second term is averaged over 
fv (j>i)fv (P2), with f V0 (Pi) =S(vi- v \ii) 5 (ui). 

The result is the modified Smoluchowski equation: 



d t c + v Q d\\c = D R d e c+ [D\\ + Ds)d\c + D ± d]_c 
1 

1 



deTex - V • C 1 • F ex 
dersp - V • C 1 • F S p , (4) 



where d\\ = u ■ V and d± = V u(u • V). The 
convective term on the left hand side of (j4]) is a trivial 
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consequence of self propulsion and describes mass flux 
along the long axis of the rod. The first three terms 
on the right hand side of the equation describe transla- 
tional diffusion longitudinal (Dm) and transverse (Dj_) 
to the rod's long axis and rotational diffusion (D B ). 
For long thin rods Dm = 2D± = D. At low density 
D = k B T a /0i and D B = 6D/£ 2 . A novel consequence 
of self-propulsion is the enhancement of longitudinal 
diffusion by Ds = w§/C||- This can be understood by 
noting that a diffusing rod performs a random walk 
with a step length x a = C~pV/3- For thermal systems the 
rod's velocity is isotropic on average and has magnitude 
Vth ~ VkisTa- In this case the anisotropy of diffusion 
arises solely from the anisotropy of the friction tensor. 
For self-propelled rods the step length along the long 
direction of the rod is enhanced, yielding an additional 
contribution to the longitudinal diffusion coefficient. 
Equivalently, longitudinal diffusion of a self-propelled 
rod can be reformulated as a persistent random walk 
where the rod has a bias ~ Vo towards steps along its long 
axis The next three terms in (@]) describe excluded 
volume effects within the mean-field approximation due 
to Onsager. The corresponding forces and torque can 
be derived from the familiar excluded volume potential 
as r ex = -deVex and F ex = -W ex , with V ex {xi) = 
k B T a c{xi,t) /^ 12 / fl2 |ui x u 2 |c(ri + £i 2 , 2 , t), with 
£12 = £1 — £2- Finally, t$p and F$p describe, within 
a mean-field approximation, the additional torque and 
force due to anisotropic linear and angular momentum 
transfer during the collision of two self-propelled rods, 

(I 8 ") = <\ (- if V (*i x u 2 )] 2 

V t sp J Jx 2 , Sl ,s 2 ,£ \ z • (£i x k) J 

X0(-Ui2 • k)S (T cont ) C(xi,t)c(x2,t), (5) 



with U12 = iii — U2- In Onsager's mean field model, two 
thin rods of length I exchange an average momentum 
(|Av|) v coU ~ k B T a ji per unit time upon collision, with 
(|Av|) ~ \Jk B T a and v coU = v th /£ ~ \Jk B T a jl. When 
rods are self propelled there are anisotropic contributions 
to both the momentum exchanged ((|Av|) ~ i>o|ui x u 2 |) 
and the collision rate (y C oll ~ t>o|ui x u 2 |/i). These yield 
the new anisotropic steric forces and torques in Eq. ([5]). 

Hydrodynamics. We now use the modified Smolu- 
chowski equation to obtain coarse-grained equations that 
describe the dynamics of the systems on wavelengths 
long compared to the length of the rods and on time 
scales long compared to the collision time. In this 
regime the dynamics is controlled by the "slow vari- 
ables" corresponding to the conserved densities (here only 
the concentration of filaments p = J~ c(x, t) ) and the 
fields associated with possible broken symmetries. In 
a liquid of self-propelled rods, both polar and nematic 
order are possible, described by a polarization vector 
P(r, t) = J~uc(x,t) and the nematic alignment tensor 
Qaflfat) = fsiuaUp - \5 a p)c(x,t), respectively. Since 
each rod has a self propulsion velocity v \i, the polariza- 
tion is also proportional to the self propulsion flow field. 
The equations for these continuum fields are obtained by 
taking the corresponding moments of the Smoluchowski 
equation (@]) and are given by 



d t p + v ^ ■V = D p SJ 2 p + D Q VV :pQ (6) 



dtP + DrP AP • Q + v V ■ Q + y Vp + A'[3(P • V)P - ^VP 2 PV P] = D b V 2 P + (D 3pl - D b )VV ■ P (7) 



P 



~ ~pin) ^ + VoF = - x " (l p ' VQ + i QV ' p+ i G + + ^ (w - \ 1)p (8) 



where F a/3 = (d a Pp + dpP a - (5 Q/3 V ■ P) and G Q/3 = 
Qa-yd-fPp + Qf3yd<yP a — SapQa-yd-yPa- All A parameters 
in Eqs. ([7]) and ((HJ) are proportional to Vq and vanish in 
the absence of self propulsion. All diffusion constants are 
enhanced by self-propulsion via additive terms propor- 
tional to Ds. Finally, we have suppressed in Eqs. (JGHSJ) 
excluded volume corrections to the diffusive terms, non- 
linear terms of second order in gradients, and corrections 
to the convective terms beyond linear in Vq. The com- 



plete hydrodynamic equations with explicit expressions 
for the various coefficients can be found in Ref. 

The stable homogeneous stationary solution of Eqs. ([6]- 
[8]) are the bulk states of the self-propelled system. Two 
such states are possible: an isotropic state, with p = 
constant, P = 0, Q a p = 0, and a nematic state, with 
p = constant, P = and Q a p ^ 0. Hard core interac- 
tions and self-propulsion modeled simply as a body force 
are not sufficient to generate a bulk polar state, with 
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FIG. 2: (color online) The polarization in a channel of width 
L for S/L = 0.2 (solid) and S/L = 0.6 (dashed). 



P ^ 0. Either shape or mass distribution asymmetry 
of the driven particles or hydrodynamic interactions are 
essential to obtain a macroscopic polar (moving) state. 
Self-propulsion has, however, a profound effect on the 
isotropic-nematic transition which occurs at the density 

2 

Pin(v ) = Pjv/(1 + 5F^t), where p N = 3/(7r^ 2 ) is the 
Onsager transition density. The transition occurs where 
the coefficient of the term linear in Q a p on the right 
hand side of Eq. changes sign, signaling the unsta- 
ble growth of nematic fluctuations. This enhancement of 
orientational order has been observed in numerical simu- 
lations of actin motor assays, where actin filaments move 
on a substrate grafted with motor proteins [§]. It arises 
from the additional torque tsp that self-propelled rods 
experience upon collision as compared to thermal rods. 
This enhances entropic ordering and aligns the rods |2ll |. 

Although no bulk polar order is possible in our sys- 
tem, self-propulsion greatly enhances the length scale 
over which polarization fluctuations decay. As a result 
boundaries pay a crucial role in self propelled systems. 
To illustrate this we consider a self-propelled 2d hard 
rod fluid confined in the channel of width L between two 
boundaries, as shown in Fig. O We assume that the 
boundaries induce polarity by forcing all rods to align 
in the same direction, i.e., P x (—L/2) = P x (L/2) = Pq. 
In this geometry the density is constant. One can eas- 
ily solve for the polarization profile across the channel 
with the result P x (y) = Po c °sh(y/<5)/ cosh(L/2S), where 
5 = s/D b /D R = £/2y/5/2 + v%/k B T is the boundary 
layer width over which the polarization penetrates in the 
channel. In the absence of self-propulsion 5 ~ £, i.e., a 
finite polarization at the boundary decays (via rotational 
diffusion) over a length scale of order I. For large self- 
propulsion velocity, 6 ~ \vo\. If L ~ S the entire channel 
is effectively polarized. We stress that numerical simu- 
lations of self-propelled rods on a substrate have indeed 
observed large correlated regions of finite polarization, 
but never an ordered bulk state. We expect that the 
boundary layer length 5 also sets the scale of correlations 
in bulk systems. Finally, as shown in [IB], Eqs. (JMHl yield 
two important properties of fluctuations in self-propelled 
systems. First, the isotropic state can support sound- 
like propagating density waves for a range of wavevec- 
tors above a critical value of Vq . Secondly, large number 



fluctuations always destabilize the homogeneous nematic 
state. We refer the reader to Ref. [l5] for a complete 
description of both results. 

In summary, we have analyzed a simple model that 
captures two crucial properties of self-propelled systems: 
the orientable shape of the particles and the self propul- 
sion. Using the tools of nonequilibrium statistical me- 
chanics we have derived a modified Smoluchowski equa- 
tion for SPP and used it to identify the microscopic origin 
of several observed or observable large scale phenomena. 

This work was supported by the NSF on grants DMR- 
0305407 and DMR-0705105. 
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